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Abstract The interference of two independent single-pho- 
ton pulses impinging on a beam splitter is analysed in a gen- 
eralised time-resolved manner. Different aspects of the phe- 
nomenon are elaborated using different representations of the 
single-photon wave packets, like the decomposition into sin- 
gle-frequency field modes or spatio-temporal modes match- 
ing the photonic wave packets. Both representations lead to 
equivalent results, and a photon-by-photon analysis reveals 
that the quantum-mechanical two-photon interference can be 
interpreted as a classical one-photon interference once a first 
photon is detected. A novel time-dependent quantum-beat ef- 
fect is predicted if the interfering photons have different fre- 
quencies. The calculation also reveals that full two-photon 
fringe visibility can be achieved under almost any circum- 
stances by applying a temporal filter to the signal. 
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1 Introduction 

Single-photon light pulses are a central ingredient of quantum- 
teleportation and quantum-communication protocols as well 
as quantum computation with linear optics Q]|2][3]|3) . Many 
of these schemes are based on the quantum interference of 
two indistinguishable photons that impinge simultaneously 
on different entrance ports of a 50:50 beam splitter. In this 
case, both photons leave the beam splitter at the same output 
port. This quantum interference effect was first observed by 
Hong, Ou and Mandel with photons from a parametric down- 
conversion source ||5]|6]E1. When the two photons had differ- 
ent frequencies, a spatial quantum-beat phenomenon was ob- 
served |8|. The two-photon quantum interference effect has 
also been observed with photons from independent down- 
conversion sources lEJl llOI . and the phenomenon has been 
used by Santori et al. II II to demonstrate the indistinguisha- 
bility of photons emitted from a quantum dot embedded in a 



microcavity. A theoretical description that applies to the in- 
terference of independent single photons has been worked out 
by Bylander et al. fHI . 

So far, all these experiments use photons that are short 
compared to the time resolution of the employed photodetec- 
tors, and temporal effects are therefore neglected. However, 
we have recently developed a novel type of single-photon 
emitter that employs a strongly coupled atom-cavity system 
ll3M14lfT5l to emit single-photon wave packets of 2 /is du- 
ration. This exceeds the typical response time of solid-state 
photodetectors by more than three orders of magnitude and 
allows to study two-photon quantum interference effects in 
a time-resolved manner. This paper therefore presents a de- 
tailed analytical treatment of the two-photon interference pro- 
cess. It takes into account the time evolution of the amplitudes 
and phases of the interfering single-photon wave packets. The 
analysis is valid for independent photons, and therefore ex- 
tends the model for frequency-entangled photon pairs from 
Steinberg et al. 1 16 1. 

While preparing this paper, we realised that apparently 
there is no commonly adopted approach to deal with free- 
running single-photon wave packets 1171 . One possibility is 
to express the quantum state of a single-photon wave packet 
as a superposition state of an infinite number of single-fre- 
quency modes, another approach uses a superposition state of 
spatio-temporal modes. To demonstrate that these approaches 
lead to the same result, we now analyse the two-photon inter- 
ference effect along both routes. 



2 Two-photon interference 

As a starting point for the following discussion, we first re- 
view the two-photon interference phenomenon in the Fock- 
state picture, i.e. we consider two single photons that impinge 
on the entrance ports labelled '1' and '2' of the beam splitter 
shown in Fig.^ This initial situation corresponds to the quan- 
tum state | li I2) = ajaJ, |0), which is obtained by applying 
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the appropriate photon-creation operators a\ on the vacuum 
state |0). The effect of the 50:50 beam splitter on the field is 
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Fig. 1 Two single-photon wave packets impinge with a relative de- 
lay of 5t on ports T and '2' of a beam splitter. The pulse duration 
is longer than the time resolution of the detectors monitoring ports 
'3' and '4' of the beam splitter. Therefore, the time difference, r, 
between two photon detections must be taken into account in the 
analysis of the photon statistics. 

usually described by the unitary transformations, 
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and 



4 = (a\ +4)/\/2 

4 = (4-4)A/2- 



If we use these relations to express the initial state in terms of 
photons created in the two output ports labelled '3' and '4', 
we immediately see that this leads to an entangled state with 
either both photons in one or the other output port, 



|iii 2 ) = 44l 



(4 + 4x4-4)io) 



i([4i 

(|2 3 4 ) 



[4nio) 

|0 3 2 4 ))/V2. 



(2) 



This simple picture qualitatively explains the observed inter- 
ference phenomena, but since the photonic modes that belong 
to the creation operators are neither specified in space, fre- 
quency nor time (albeit they are assumed to be identical for 
both photons), quantitative results for photon wave packets 
cannot be obtained from this approach. 



2.1 Space-time domain 

To obtain more insight, we now consider a description of 
the arriving single-photon wave packets in the space-time do- 
main. The free-running photons are not subject to boundary 
conditions that restrict the photon creation and annihilation 
operators, a£ and a k , to specific frequency modes. Therefore 
we are free to define these operators in such a way that they 
create or annihilate photons in arbitrarily chosen modes (la- 
belled by the index k) in space and time, which we define by 
the one-dimensional spatio-temporal mode functions 



U^t) = e k {t-z/c)e- l ^ z ' c \ 



(3) 



Beside an arbitrary phase evolution, </)/ c (t), the mode func- 
tions incorporate amplitude envelopes, £&(£), which are as- 
sumed to be normalised so that J dt\ek(t)\ 2 = 1. By placing 
the beam splitter at z = 0, we can omit the spatial coordinate, 
z, in the following. Now we assume that the whole Hilbert 



space is spanned by an orthonormal set of spatio-temporal 
modes, (k(t), so that the electric field operators read 

E + (t) = ]TCfc(<H and £~(t) = £tf(*)4- 

k k 

(4) 

Provided a photon is present in mode i, the probability to 
measure it at a given time t then is 

Pi(t) = (uiE-m+m,) = a(t)Q(t) = hm 2 , (5) 

which is solely defined by the amplitude envelope of the pho- 
ton, ej(t). To analyse the effect of the beam splitter, we now 
attribute field operators, E^ 2 3 4 , to each of the four input/out- 
put ports. For the two input ports we only consider the occu- 
pied modes which are described by the mode functions £i (t) 
and (,2{t), respectively. This restricts the Hilbert subspaces 
of the input ports to single modes. Therefore the field op- 
erators that belong to the two input ports can be written as 
= Ci{t) a i an d E% = C2{t) a 2, respectively. Using these 
operators, the effect of the beam splitter is now described as 
a transformation between input and output modes: 



E+(t) = [E+(t)+E+(t)]/V2 
= [Ci(t)ai + C 2 (f)a2]/V2 

Ei(t) = [E+(t)-E+(t)]/V2 
= [CiWai-C 2 (t)a 2 ]/v / 2 



(6) 



Note that this is a more general form than the simple opera- 
tor relation Eq. Q, which is equivalent to Eq. (|6j if all mode 
functions, are equal. 

2.7.7 Joint photon-detection probability: To obtain the prob- 
ability for photon detections in the output ports '3' and '4' at 
times to and to+r, respectively, we have to apply the field op- 
erators Ef- (to ) and Ef (to + r). With two photons impinging 
on the two input ports in modes £1 (t) and £2 (t), the incoming 
state is |^i„) = 44|0), an d the joint probability for photon 
detections reads 

Pjoint{tOi T ) — 934^0, £0 + T) = 

(0\ ai a 2 E^(t Q )E^(t Q + T)E+(to + T)E+(t )a\4\0). 

(7) 

With the field-operator relations defined in Eq. @ and the re- 
striction to single-photon states, i.e. using (l|aW = aa\l) = 
0, one easily verifies that Eq. Q leads to 



P 



joint (^0 i 



T\ = j|Cl(*0" 



-T)C2(io)-C2(io + T)Ci(*o)| 2 - (8) 



The most striking feature of this joint probability function is 
that it vanishes for r = 0, no matter how different the mode 
functions of the two incoming photons are. Moreover, if the 
phase relation between £1 and (2 is lost for r 3> t c , where 
t c is the mutual coherence time of the incoming photons, the 
interference term in Eq. (JHJl, averaged over an ensemble of 
different photon pairs, is zero, and Pjdnt reduces to 



P, 



joint 



(*o,r) 



r»T c 



1 



(Pl{t )P2(t 



P2(t )Pl(to + T)) 

(9) 
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2.1.2 Photon-by-photon analysis: Despite the fact that the 
full interference phenomenon is described by Eq. l|8}, it is 
hard to acquire an intuitive understanding of the interference 
process from the two-photon approach. This is not the case 
if one analyses the effect step-by-step. We therefore use the 
same spatio-temporal mode functions as before and ask now 
for the probability to detect a photon in output port '4' at time 
to + t conditioned on a photon detection in output port '3' at 
time to. The probability to find a photon in port '3' reads 

iM«0) = (^n\E^{t )E+(t )\1r m ) = 

(l 1 l 2 |(4-(*o)+ ^(io))|(A + (io) + ^ 2 + (io))|lil2> - 
i(|ex(*o)| 2 + |e 2 (* )| 2 ). 

(10) 

Note that this expression includes no interference term, since 
the initial product state is composed of single-photon Fock 
states that have no relative phase. The state \&cond), condi- 
tioned on the detection of a photon in port '3' at time to, is 
obtained by applying the field operator £3 (to) to the incom- 
ing state, 1 1 1 12) - After re-normalisation, this leads to 



4 + (*o)|lil; 



C2(*o)|li0 2 ) + Ci(*o)|Oil 2 ) 



id)- 



(11) 

The conditioned state, \W con d), describes a single photon that 
either impinges on port '1' or port '2' of the beam splitter, 
with the amplitude and phase relations between these two 
modes given by C^ito) and (i(*o)- This photon now gives 
rise to classical interference fringes that start in phase at time 
to, i.e. the second photon is found in the same port as the 
first photon if it is detected at the same instant. However, the 
phase evolves in time according to £1 (t) and (2 (t). Therefore 
the probabilities to detect the second photon in output ports 
'3' or '4' at time to + t read 

^3,4^0 + t) = {^amd\E^ A {to + T ) E 3 A {t + T)\<P cond ) 

_ |Ci(*o + r)Ca(to) ± C 2 (io + r)Ci(t )| 2 



2(|ei(to)| 2 + h(*o)| 2 ) 



(12) 



Together with Eq. dl0> . the joint photon-detection probability 
for a first detection in output port '3' at time to and a subse- 
quent photon detection in output port '4' at time to + r is 

Pjoint(t ,T) = P 3 (to)Pi(to + t) 

= ilCi(*o +r)C 2 (<o) -C 2 (to + r)Ci(to)| 2 . 

(13) 

Result Jl 3I > is identical to expression (|8j. However, the step- 
by-step analysis reveals the equivalence between classical one- 
photon interference and two-photon interference conditioned 
on the detection of a first photon. Note that the field operators 
at the output-ports commute, i.e. [E^ (to), E± (to + t)] = 0. 
Therefore the sequence of operators is irrelevant, and for neg- 
ative detection-time delay, r, which corresponds to a reverse 
order of photon detections, the same results are obtained. 



2.2 Single-frequency field modes 

A standard way to attack the problem of two interfering pho- 
tons is based on the decomposition of single-photon wave 
packets into an infinite number of single-frequency field modes. 
To do so, a single photon in one of the two spatio-temporal 
input modes (i = 1, 2) is expressed as one quantum of ex- 
citation occupying a superposition of all possible frequency 
modes, 

|1<) = [ dco $i(w) at(w) |0), (14) 



where the spectral amplitude of the single-photon pulse, <Pi(uj), 
is assumed to be normalised so that J dcu $* (u})4>i(u>) = 1. 
The continuum of spectral fields allows one to express the 
time-dependent electric field operators as 



Et(t) 



1 



2tt 
1 



dio 



t a,i(uj) and 



27T 



4)e <ui aj(w). 



(15) 



To evaluate the effect of the field operator, Ef(t), acting on 
the single-photon state |lj) defined in Eq. (I14t . we make use 
of the fact that the Fourier transformation of the spectrum, 
<Pi(io), equals the spatio-temporal mode-function, 



dco 'Pi (u>) e 



(16) 



which leads to 



E+(t)\W = TkN dw ^ e'^H *i(w)aj(d>) |0) 



o) = Ci(t)|o). 

(17) 

Therefore, one easily verifies that the probability to find the 
photon at time t reads 



Pi(t) = (h\Er(t)E+(t)\U) 

= {0\Q(t)Q(t)\0) = Q(t)Ci(t), 



(18) 



which is equivalent to Eq. (|5jl, since the two approaches were 
mapped onto one another (see Eq. dl6ll). 

In the field-mode picture, the effect of the beam splitter 
is described by the transformation relations in Eq. Q, and 
the quantum state of the two incoming photons, 1 1 1 1 2 ) , is ex- 
pressed as a product state of field-mode superposition states, 
as they are defined in Eq. dl4> . Therefore we make use of 
the relations in Eq. ([3 and replace a}(wi) and 02(^2) by the 
respective sums and differences of a 3 (uj) and a\(ui), which 
leads to 

II1I2) = JJ duJ\duJ2 $i(wi)$ 2 (£«>2) ai(fc>i)4(cj 2 )|0) 
= jj du>idu> 2 ^1(^1)^2(^2) 

x±[4(wi) + <4(wi)][4(w 2 ) - 4(^2)]|o). 

(19) 

We now use Eq. to calculate the joint detection-probability, 



Pjoint{to, t) — 

(hl 2 \Ez(t )E7(t + r)E+(t + r)4+(t )|l 2 l 2 ) 



(20) 
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for two photons detected in output ports '3' and '4' at times 
to and to + r, respectively. We use Eq. (1191 to express the 
state | I1I2) in terms of the creation operators a\ and a\ to 
obtain the right-hand-side of the operator-bracket in Eq. (I20> . 

Ei{t + T)E+{t )\l x lz) = 

X [6(u>3 — 0J 2 )S(ui4~ UJ\) — S(uj 3 — ^1)6(014 — UJ 2 )]\0). 

(21) 

We now make use of the Fourier-transformation relation be- 
tween $i(ui) and d(t) from Eq. (fT6l to evaluate the above 
integral, which leads to 

E+ (t + t)E+ (t )|lil2> = 

1 (22) 

2[Ci(* + r)Ca(*)-Ci(*)C 3 (t + r)]|0), 

so that the joint photon-detection probability finally reads 

Pjoint(t0,T) = i|Cl(*0 + 7-)Ca(*0) " Cl(*0)C2(*0 + r)\ 2 , 

(23) 

which is again equivalent to Eq. (|8}. Therefore we conclude 
that the different approaches of modelling single-photon wave 
packets all lead to the same expression for the joint two- 
photon detection probability, Pj int(to, T )- But note that the 
standard way of using single-frequency field modes is not well 
adapted to the problem and requires a lengthy and cumber- 
some calculation. This is not the case in the space-time do- 
main, where the joint detection probability can be calculated 
in a fast and elegant way (section I2.1.H . However, physical 
insight is only obtained from the photon-by-photon analysis 
(section 12.1.2k which reveals that the two-photon quantum 
interference is in fact equivalent to a classical interference 
phenomenon once a first photon detection has taken place. 
In the collapsed state, the relative phases and amplitudes of 
the two interfering modes are fixed. Therefore an oscillating 
quantum-beat signal is expected if the two spatio-temporal 
mode functions, &(<;) an d evolve at different rates. 

3 Quantum beat 

To illustrate the quantum-beat signal in the interference of 
two photons with slightly different frequency, we now con- 
sider the case of two equally long Fourier-limited Gaussian 
single-photon pulses that impinge on a 50:50 beam splitter 
with a relative delay time St and carrier-frequency difference 
A = uj 2 — uj\. These pulses are described by the two spatio- 
temporal mode functions 

&(t) = </27^exp(-(t - St/2) 2 - i(w - A/2)t) 
and (24) 
( 2 (t) = ^/2/Vexp(-(i + St/2) 2 - i(u + A/2)t), 

with uj = (u>i + uj 2 )/2. For the sake of simplicity, the time, 
t, the photon delay, St, and the detection-time delay, r, are 
expressed in units of the pulse duration, St (half width at 1 /e 




Fig. 2 Joint probability, P^hv, for the detection of photon pairs as 
a function of the relative delay between the two impinging photon 
wave packets, St, and the delay between photon detections, r for the 
three frequency differences A — 0; n/2; 3ir. All times and frequen- 
cies are normalised by the photon-pulse duration, St, and, hence, are 
dimensionless. 



maximum), while the frequencies, lo and A, are expressed in 
units of 1 /St. 

With the spatio-temporal modes defined in Eq. i24\ . a straight- 
forward evaluation of Eq. © leads to the joint photon-detection 



Time-Resolved Two-Photon Quantum Interference 



5 



probability 

Pjotnt(to,T, St, A) = 

cos1i(2t<5t) — cos(tA) , , , \ r , „ 
* '- * '- x cxp(-4i (to + t) - St 2 - 2r 2 ), 

7T 

(25) 

which depends on the photon detection times, to and to + r. 
To obtain the probability of detecting two photons in the ports 
'3' and '4' with a time difference r, we integrate over all 
possible values of to, which leads to 



P 2Hv {t,5t, A) = gi%T)\ ST ,A 

J dt Pjoint(to,T,8T, A) = 

cos1i(2t(5t) — cos(tZI) 



(26) 



2^ 



x exp(-<5r 2 - t 1 ). 



This joint two-photon detection probability is displayed in 
Fig.|2]as a function of the photon-delay, St, and the detection- 
time delay, r, for three different values of the frequency dif- 
ference, A. 

The case with A = is depicted in Fig.|2ja). As ex- 
pected, no coincidences are found if the photons arrive si- 
multaneously, i.e. for St = 0. This situation is well explained 
by the Fock-state model of two interfering indistinguishable 
photons. Moreover, even if the two photons are delayed with 
respect to each other, we see that simultaneous photon de- 
tections (with t — 0) never occur. This is well explained by 
Eq. which yields that F 'joint (to,T — 0) = 0, no matter 
how different the mode functions are. In particular, the sec- 
ond photon leaves the beam splitter through the same port 
as the first one, provided no dephasing takes place, i.e. for 
|Ci(*o + r)Ca(*o) - Ci(MC 2 (<o + t)\ « 1. 

Figure|2jb and c) show situations where the two imping- 
ing photon wave packets have a frequency difference of A = 
ir/2 or A — 3ir, respectively. Most remarkable is that simul- 
taneous detections (with r = 0) still do not occur, although 
the two photons are now distinguishable. Moreover, for si- 
multaneously arriving photon wave packets (with St = 0), 
the coincidence probability now oscillates as a function of t 
with the frequency difference A. We emphasise that this beat 
note always starts at zero for t = 0, due to the in-phase start- 
ing condition imposed by the detection of the first photon. 

It follows that the coincidence probability remains zero 
even in case of an inhomogeneous broadening of the fre- 
quency differences, A, To illustrate this, we assume a Gaus- 
sian frequency distribution, 



1 



Slu^/tt 



exp(-(/y<W) 2 ). 



(27) 



where Slo defines the bandwidth of the spectrum (half width 
at 1/e maximum). The integral of this frequency distribution 
is normalised to one, so that we can use it to calculate the pho- 
ton detection probability for simultaneously arriving photon 



0.25 



0.2 



"S 0.15 



0.1 



■a cos 



o 
U 



- envelope (5(0 

- 5co = 3k [Sr 1 

- ■5co = 7i [5t" 

- 5tfl = 




-2-10 1 2 

Detection-time difference x [5t] 

Fig. 3 Time-resolved two-photon detection probability, Pi n h, as a 
function of the detection-time difference, r, in case of simultane- 
ously arriving equally long single-photon wave packets (St = 0) 
for different inhomogeneous spectral widths, 8u>. The solid line rep- 
resents non-interfering photons and therefore reflects the amplitude 
envelope of the photonic wave packets. 




3jc >^ 



Fig. 4 Total probability for the detection of photon pairs as a func- 
tion of the relative delay between the two impinging photon wave 
packets, St, and the inhomogeneous spectral width, Slo, of the inter- 
fering photons. 

wave packets (with St = 0) by evaluating 

P mh {T,Suo) = jdA f Sa {A)P 2hv (T,ST = 0,A) 



exp(— r 2 ) 
20? 



1 — exp — 



\2/Su 



(28) 

Figure [5] shows that the inhomogeneous broadening of the 
photonic spectrum leads to a dip in the two-photon detection 
probability with width 2/Su> (half width at 1/e dip-depth). We 
emphasise that this dip in Pi n h reaches zero. Therefore it is 
possible to achieve perfect two-photon coalescence if a tem- 
poral filter is applied that restricts the detection-time differ- 
ence to |t| <§; 2/ Slo. 

To connect our results to previous experiments with ultra- 
short photons, we now calculate the total two-photon coinci- 
dence probability by integrating over the detection-time de- 
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lay, r. In case of inhomogeneous broadening, as discussed 
above, this leads to 

P tota i{6T,6tu)= [[ dTdAf SuJ {A)P 2hu {r,5T,A) 



2 V4 + Su 2 ' 

The total two-photon coincidence probability is shown in Fig.|4] 
as a function of the arrival-time delay, 5t, and the inhomoge- 
neous spectral width, Sui. In this case, it is evident that the 
two-photon coincidence dip reaches zero only if Sco = 0, i.e. 
only if the two photons are mutually coherent. Increasing the 
bandwidth has no effect on the width of the dip, but only de- 
creases its depth. Therefore spectral filtering is mandatory if 
the time resolution does not allow to filter photon-detection 
times. 
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4 Conclusion 

A rich substructure is expected in two-photon interference 
experiments if the photons can be detected with high time 
resolution. For example, distinguishable photons of different 
frequencies are expected to give rise to a pronounced oscilla- 
tion of the joint photon-detection probability behind the beam 
splitter. The ability of generating long single-photon pulses 
by means of a strongly-coupled atom-cavity system opens up 
the possibility to 'zoom' into the quantum state of a single- 
photon pulse by measuring two-photon quantum interference 
fringes as a function of time. 
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